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ABSTRACT 

McKay et al. (2002) have recently used measurements of weak galaxy-galaxy lensing 
in the Sloan Digital Sky Survey to estimate the cross correlation between galaxies 
and the projected dark matter density field. They derive a relation between aperture 
mass within a radius of 260ft.~ 1 kpc, M260, and lens galaxy luminosity, that does not 
depend on galaxy luminosity, type or environment. In this paper, we study the cross- 
correlation between galaxies and dark matter using galaxy catalogs constructed from 
a high-resolution N-body simulation of a ACDM Universe. We show that our simula- 
tions reproduce the McKay et al. results reasonably well. In the simulation, M260 is 
approximately equal to the halo virial mass for L* galaxies. Mim overestimates the 
virial mass for fainter galaxies and underestimates it for brighter galaxies. We use 
the simulations to show that under certain circumstances the halo virial mass may 
be recovered by fitting an NFW profile to the projected galaxy-mass correlation func- 
tion. If we apply our method to the observations, we find that L* galaxies typically 
reside in halos of ~ 2 x 1O 12 /i -1 M , consistent with halo masses estimated from the 
observed Tully-Fisher relation. In the simulations, the halo virial mass scales with 
galaxy luminosity as L for central galaxies in halos and for galaxies in low-density 
regions. For satellite galaxies, and for galaxies in high-density regions, there is no 
simple relation between galaxy luminosity and halo mass and care must be exercised 
when interpreting the lensing results. 

Key words: dark matter - gravitational lensing - large-scale structure of the universe 
- galaxies: halos - methods: statistical 



1 INTRODUCTION 

In a recent study, McKay et al. (2002, herafter M2002) used 
the signal of galaxy-galaxy lensing detected in the Sloan Dig- 
ital Sky Survey (SDSS) to infer the dark matter mass dis- 
tribution around bright galaxies. What McKay et al. mea- 
sured from the data was the shear field 7+, inferred from 
the distorted images of background galaxies around lensing 
galaxies. From the theory of weak lensing (e.g. Schneider 
et al., 1999), the shear field is related to the surface mass 
density contrast AE+ as: 



7+(^)S c 



E(< R) - Z(R) = AE 4 



(1) 



where E(< R) is the mean surface density within radius R, 
E(i?) is the mean surface density at radius R, and E cr it is 
the critical density determined by the geometry of the lens- 
source system. Defined this way, AE+ is a measure of the 
cross correlation between the lens galaxies and the projected 



mass density. We will therefore refer AE4 
mass correlation function (GMCF). 



as the galaxy- 



The GMCF depends on the matter distribution around 
lens galaxies. The observed GMCF can therefore be used to 
infer the extent of dark haloes around galaxies (e.g. Hudson 
et al. 1998; Fischer et al. 2000; Smith et al. 2001; Wilson et 
al. 2001; Guzik & Seljak 2002). Using the SDSS data, McKay 
et al. (M2002) found that: (i) the shear field 7+ could be 
measured with reasonable signal-to-noise out to a projected 
radius ~ l/i~ 1 Mpc; (ii) AE+(J?) could be fit by a power 
law, with the singular isothermal sphere an acceptable model 
for the distribution of dark matter around galaxies; (iii) the 
mass M260 (defined as the mass projected within an aperture 
of radius 260 /i _1 kpc) inferred from the observed GMCF was 
strongly correlated with the luminosity of lens galaxies and 
in the redder photometric pass-bands, the galaxy luminosity 
was proportional to M260, with a proportionality constant 
of order 100; (iv) the ratio of aperture mass to luminosity in 
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the red photometric bands did not depend significantly on 
galaxy luminosity, type and environment. 

Because current models of galaxy formation based on 
the cold dark matter (CDM) cosmogony predict the exis- 
tence of extended dark haloes around galaxies, the results 
of M2002 provide important constraints on models of galaxy 
formation. 

The observed GMCF cannot be interpreted simply in 
terms of individual halo profiles. There are three other fac- 
tors that can affect it. First, the observed GMCF is an aver- 
age of mass density fields around different galaxies. Second, 
the GMCF depends on how lens galaxies are distributed 
within dark haloes. Third, more than one halo can con- 
tribute to the projected density field around a given galaxy. 
These effects must be quantified in order to use the observed 
GMCF to infer the properties of galaxy haloes. 

In this paper, we use a galaxy catalogue constructed 
from a high-resolution N-body simulation to understand 
the observed galaxy-galaxy lensing results. We first examine 
whether or not the simulation results agree with the obser- 
vations. We then illustrate how the three effects mentioned 
above influence the GMCF. Finally, we discuss what kind of 
information about dark haloes it is possible to infer from the 
GMCF. The paper is organized as follows. In Section 2 we 
describe the simulation data. Calculations of the GMCFs 
for various samples are presented in Section 3. In Section 
4, we investigate the relation between the masses derived 
from the GMCFs and the actual masses of dark haloes in 
the simulation. We then use these relations to estimate the 
halo masses of observed lens galaxies. Finally, in Section 5, 
we summarize and discuss our results. 



2 THE SIMULATIONS 

In this paper, we use the results of the GIF simulation of 
the ACDM model [with (fi , A, T, cr 8 ) = (0.3, 0.7, 0.21, 0.90)] 
carried out by Kauffmann et al. (1999) using codes from the 
Virgo Consortium (Jenkins et al., 1998; the data are publicly 
available 0.) The simulation evolves 256 3 cold dark matter 
particles in a periodic cube of side length 141.3h _1 Mpc, with 
particle mass 1.4 x 10 10 hT 1 Mq, and with a force softening 
rj = 20/i _1 kpc. Dark haloes were selected from the simu- 
lation using a "friends-of-friends" (FOF) algorithm with a 
linking length of 0.2 times the mean inter-particle separa- 
tion, and galaxy catalogues were constructed by populat- 
ing dark haloes with galaxies according to a semi-analytic 
model of galaxy formation described in Kauffmann et al. 
(1999a). Model galaxies are assigned luminosities and mor- 
phological types according to their star formation histories 
and bulge-to-disk ratios. Previous investigations with these 
simulated catalogues demonstrated that the model is quite 
successful in matching the observed properties of the galaxy 
population, both at the present day and at high redshifts 
(Kauffmann et al., 1999a, b; Diaferio et al., 1999; 2001). 

For this study, we use new publically available galaxy 
catalogues constructed from the GIF simulations giving pre- 
dictions for the SDSS photometric bands. We present results 
only for the z-band. Calculations were also made for the r 

* http:/ /www. mpa-garching.mpg.de/Virgo/ 




z' band absolute magnitude 



Figure 1. The distribution of GIF galaxies with respect to the 
z-band absolute magnitudes. Results are shown separately for 
central galaxies and all galaxies. 



and i-bands, but the results are very similar. Fig. 1 shows 
the luminosity distribution of the model galaxies. Results are 
shown for 'central galaxies' (which are the brightest galax- 
ies located at the centers of dark matter halos) and for 'all 
galaxies' (which include both central galaxies and satellites). 
As can be seen, the bright end of the luminosity function is 
dominated by central galaxies, while satellite galaxies dom- 
inate the faint end. Because of the mass resolution of the 
simulation, predictions for galaxies fainter than —19.5 mag 
are uncertain. Such galaxies are excluded from our analyses. 

3 ANALYSES 
3.1 The GMCF 

As explained in Section 1, the observed shear fields around 
galaxies can be used to derive the cross-correlation between 
lens galaxies and the projected mass density field. In this 
paper, we calculate the GMCF directly from the mass and 
galaxy distributions in the simulations. We thus neglect all 
the uncertainties involved in reconstructing the mass-density 
field from the weak lensing observations. 

We first project the positions of galaxies and dark mat- 
ter particles onto a plane (chosen to be one surface of the 
simulation box) . For each galaxy, we estimate the mean sur- 
face density contrast of dark matter within rings of different 
radii around the galaxy. The GMCF as a function of radius, 
AE+(i?), is obtained by averaging over a given sample of 
galaxies. Since by definition the background surface density 
is subtracted from AE+(i?), the derived GMCF is indepen- 
dent of the depth of the projection (here equal to the sim- 
ulation box size), so long as this depth is much larger than 
the correlation length of dark matter. 

The GMCF is estimated for two different kinds of lens 
galaxy samples: 
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• AG-sample, which contains all simulated galaxies in 
some luminosity range; 

• CG-sample, which contains only central galaxies in 
some luminosity range. 

To quantify the GMCF, we fit it by simple models. Fol- 
lowing M2002, we first fit a GMCF with a power-law, 



AE+ = A {R/lh- 1 Mpc) b MVl pc~ 



(2) 



where A and b are the fitting parameters. Motivated by the 
fact that the density profiles of CDM haloes are well de- 
scribed by the NFW profile (Navarro, Frenk & White 1997) , 
we also fit the GMCF by assuming the halo profile to have 
the form: 



p{r) 



8p 



(r/r s )(l + r/r s y 



(3) 



where 5 is a dimensionless density amplitude, r s is a charac- 
teristic radius, and p is the mean background density. This 
profile is characterized by a virial radius V2oo (defined to be 
the radius within which the average mass density is equal to 
200p), and a dimensionless concentration parameter c (or a 
scale radius r s = r2oo / c) . The amplitude 8 is related to c by 



200C 1 



3[m(l + c) + c/(l + c) 



(4) 



The properties of the projected NFW profile are given in 
Appendix A. 



We fit the GMCFs by equation (A6), with c and r2oo 
as free parameters. Note that we are not fitting a single 
dark halo with the NFW profile (the GMCF is a superpo- 
sition of different haloes), and so the values of c and r2oo 
are not related to the concentration and radius of a single 
halo but the mean of many haloes. For both the power-law 
and NFW model, we restrict the fits to the range of radii 
between 0.05/i _1 Mpc and lft _1 Mpc. 

In Fig. 2 we show the results of the GMCF for both AG 
and CG samples brighter than — 19. 5 m . These results are ob- 
tained using all simulated galaxies down to this magnitude 
limit and do not take into account the selection function 
in the observed data. Observational selection effects will be 
considered in subsection 3.3. As one can see from Fig. 2, the 
GMCFs for the AG and CG samples have rather different 
behaviour. For the CG sample, the GMCF is strongly peaked 
at small radii and tends towards zero at large radii. In this 
case, the GMCF profile is well-described by the simple mod- 
els based on the projection of halo profiles. Not surprisingly, 
the NFW model gives a better fit than the power-law model. 
In contrast, the GMCF for the AG sample is less centrally 
peaked and tends towards a constant non-zero value at large 
radii. This occurs because many satellite galaxies are lo- 
cated at the outskirts of massive halos in the simulations. 
As demonstrated in Appendix B, the reduction in the am- 
plitude of the GMCF on small scales is due to the separa- 
tion between satellite galaxies and halo centers. The posi- 
tive value at large radii is due to the extended distribution 
of mass in massive haloes. Note that the GMCF profile for 
the AG sample does not follow the simple models based the 
projections of halo density profiles. 
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Figure 2. The GMCFs for the CG (squares in the upper panel) 
and AG (squares in the lower panel) samples constructed from the 
GIF galaxies with z'-magnitudes brighter than — 19. 5 m . The solid 
curve is the fit to the NFW model, with the fitting parameters 
c and r2oo (in h~ 1 Mpc) given in N(c,r2oo). The dotted curve is 
the fit to the power-law model, with the fitting parameters A and 
b given in P(A, b). 



3.2 Dependence on galaxy properties 

Fig. 3 and 4 show how the predicted GMCFs vary as a func- 
tion of galaxy absolute magnitude. For bright magnitudes, 
most galaxies are central galaxies, and so both AG and CG 
samples give similar results. For galaxies fainter than about 
— 21 m the AG and CG samples behave differently, because 
the contribution of satellite galaxies becomes important in 
the AG sample (for example, over one third of the galax- 
ies with magnitudes in the range — 21 m - — 20 m are satel- 
lite galaxies, see Fig. 1). There is a clear trend that more 
luminous galaxies exhibit a stronger cross-correlation with 
the mass. This is expected, because luminous galaxies reside 
preferentially in more massive haloes. As we will show later, 
this result is reproduced quantitatively in the observational 
data. 

As one can see, the GMCFs are well described by the 
NFW model for for central galaxies and for bright galax- 
ies. The fit to the power-law model is not as good, but the 
difference is not big enough to be detectable in the current 
SDSS data. Note that for AG samples fainter than -21 m , 
the contamination by satellite galaxies becomes so impor- 
tant that no simple model can fit the GMCFs (see Fig. 4). 
These results show that the interpretation of the GMCF is 
much simpler for bright central galaxies. 

In our simulation, the morphological classification of 
galaxies is made according to their B-band disk-to-bulge ra- 
tios. If Ms^uige — Ms,totai < 1 mag, the galaxy is classi- 
fied as an early-type (elliptical or SO) galaxy; otherwise it 
is classified as a late-type (spiral) galaxy (see Kauffmann et 
al., 1999a for details). Due to the limited simulation resolu- 
tion, morphological types are only assigned to galaxies with 
haloes containing more than ~ 100 particles. About 7850 
(13500) galaxies in the CG (AG) sample are classified as el- 
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Figure 3. The GMCFs and the model fits for 2-band galaxies in 
different luminosity ranges (CG samples). The best fit parameters 
are given as N(c, r2oo) f° r the NFW model and as P(A, b) for the 
power-law model. The errorbars are the scatter of the GMCFs 
obtained from three projections of the simulated density field. 



Figure 5. The GMCFs for elliptical (squares) and spiral (circles) 
galaxies. The upper panel is for CG samples, while the lower panel 
is for AG samples. Curves show the fits to the NFW model. 
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Figure 4. The same as Fig. 3, but for AG samples. For galaxies 
fainter than — 21 m , both power-law model and the NFW model 
fail to give an acceptable fit. 



liptical galaxies, while about 53800 (80000) galaxies in the 
CG (AG) sample are classified as spiral galaxies. In Fig. ^ we 
compare the GMCFs of elliptical and spiral galaxies. Once 
again results are shown for all simulated galaxies without 
taking into account the observational selection function. As 
can be seen, the cross-corrrelation between early-type galax- 
ies and dark matter is stronger than for late-type galaxies. In 
the simulation, these results come from the fact that early- 
type galaxies are preferentially located in massive haloes. 



3.3 Incorporating observational selection effects 

In order to compare model predictions with observational re- 
sults, we should select galaxy samples in the same way as in 
the observations. The galaxy samples analyzed in M2002 are 
apparent - magnitude limited, and so fainter galaxies have 
smaller probability of being included. In this section, we se- 
lect model galaxies according to the observed luminosity dis- 
tribution of the lens galaxies (see figure 2 in M2002) . There 
are some subtle effects associated with a magnitude - limited 
sample. For example, lens galaxies with different luminosi- 
ties have different redshift distributions, and so their lensing 
strengths (i.e. the value of £~j t ) are different. Furthermore, 
the angular radius corresponding to a given projected phys- 
ical radius is larger at a smaller distance, and so the number 
of lensed sources (N BOUIce ) within a given projected physical 
radius is larger for closer galaxies. In the analyses of M2002, 
lens galaxxies are weighted by Afsources * S~? t . To mimic this 
weighting, we use the curve of £~; t given in M2002 and as- 
sume N source oc r~ 2 , where r is the angular-size distance 
to a lens. We calculate a mean weight for galaxies with a 
given luminosity by averaging N aources * £~; t over the vol- 
ume within the distance limit which corresponds to the lu- 
minosity in consideration, assuming that the distribution of 
such galaxies within the volume is uniform. We then assign 
to each model galaxy a mean weight corresponding to its 
luminosity. Note that the iVsource used in M2002 takes into 
account other selection effects, such as the exclusion of lens 
galaxies around which the distributions of source galaxies 
are too asymmetric (see M2002 for details). Such selections 
may change the N aource oc r~ 2 relation but are not taken 
into account in our model. 

Fig.^ shows the results of the GMCF for a sample of 
galaxies selected in the z-band. It turns out that most of the 
galaxies in the selected samples are luminous central galaxies 
and as a result, the GMCF can now be fitted reasonably well 
by the NFW model. Compared with the GMCFs derived 
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Figure 6. The GMCF and the model fit for the sample con- 
structed with the SDSS selection function. The solid line is the 
fit to the NFW model, while the dashed line is half of the model 
fit. For comparison, the GMCFs derived by M2002 from the com- 
bined g'~ r'- and i'-band images are shown as squares with error- 
bars. Note that M2002 only used the images in these three bands 
to measure the galaxy shapes. Although the predicted GMCF has 
the same shape as the observation, the predicted amplitude is too 
high by a factor of ~ 2. 

by M2002, our predicted GMCF has a similar shape but a 
higher amplitude which corresponds to more massive halos. 

When observational selection is applied to the GIF sim- 
ulation, about 3200 galaxies are classified as elliptical galax- 
ies, while about 7500 galaxies are classified as spiral galax- 
ies. The fraction of the early-type galaxies in the simula- 
tion is lower than that in the SDSS sample used by M2002 
(E : S ~ 1 : 1). This may be due to the fact that the defi- 
nitions used to define early- and late-type galaxies are dif- 
ferent in the GIF simulation and in the observations. As 
mentioned above, galaxies in the GIF simulation are as- 
signed morphological types according to their bulge/disk ra- 
tios, while galaxies used in M2002 are classified according to 
their color, concentration and asymmetry. In Fig.j?] we com- 
pare the GMCFs of elliptical and spiral galaxies. Early-type 
galaxies have systematically higher GMCF than late-type 
galaxies, consistent with the observational results of M2002. 
The amplitude of GMCF is higher than the observational 
results. We will discuss this in more detail in the next sec- 
tion. 

We have also studied the effect of environment on the 
GMCF. We count the number of neighbours with projected 
distances less than l/i _1 Mpc from each galaxy and divide 
galaxies into two high- and low-density subsamples contain- 
ing equal numbers of objects. This procedure is is similar 
in spirit, but somewhat simpler than the one adopted by 
M2002, who defined the local density at the position of a 
galaxy by the inverse of the area of the polygon obtained 
from a Voronoi tesselation of the galaxy distribution on the 
sky. 

The mean z'-band luminosities for the high and low- 
density samples are 3.3 and 2.0 xl0 10 h" 2 L Q . The GMCFs 
for these samples are shown in Fig.B. Galaxies in low-density 



Figure 7. The GMCFs for elliptical (squares) and spiral (circles) 
samples with the SDSS selection function. Curves are fits to the 
NFW model. 
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Figure 8. The GMCF for galaxies in high- (squares) and low- 
density (circles) regions, selected with the SDSS selection func- 
tion. Curves are fits to the NFW model. 

environments are mainly central galaxies and as a result, 
their GMCF is fit very well by an NFW halo. Galaxies 
in high-density environments occur in groups and clusters. 
Many of them are not central galaxies- this is why the NFW 
model is a considerably worse fit. 



4 THE MASS-LUMINOSITY RELATIONS 

One important goal of studying galaxy-galaxy lensing is to 
derive the relation between luminosities of galaxies and the 
masses of their host haloes. In this section, we study whether 
or not the observed GMCF can be used to infer such a mass- 
luminosity relation. 
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Figure 9. The M260 _ l um i n osity relation for galaxies in the GIF 
simulation. Open circles are based on the fit to the NFW model, 
while solid circles are based on the fit to the SIS model. The 
squares with errorbars are the results from M2002. 



4.1 The M 26 o-L Relation 

In their paper, M2002 attempted to estimate the halo mass 
within a fixed aperture of radius R = 260/i~ 1 kpc. The rea- 
son for this choice is that the GMCF at a much larger radius 
is likely to be affected by projection effects, and at small ra- 
dius, the GMCF is difficult to measure observationally (see 
M2002 for detailed discussion). 

In this section we follow the procedure adopted by 
M2002 and fit the GMCF to a projected singular isothermal 
sphere (SIS) model to obtain M260 for the galaxy samples 
we have selected out of the simulations: 



A + X(R)sis = 



SIS 



2GR 



(•») 



where a v is the line-of-sight velocity dispersion. We have also 
fit the GMCF to a NFW profile, which gives the parameters 
r2oo and c (see Subsection 2.2). The aperture mass M260 is 
then the integral of AS+ within a radius 260/i -1 kpc. 

Fig.^ shows M26O as a function of the mean luminosity 
of galaxies. Both the SIS and NFW models give very similar 
aperture masses. For comparison, we also display the M260- 
L relation given by M2002. The predicted trend of mass with 
luminosity is similar to the observation, but the predicted 
M260 is about a factor of 2 larger than the observed value. 

In Figjl^ we show the predicted M260 / L ratios for var- 
ious simulated samples, and compare them with the mea- 
surements of M2002. Note that the observational selection 
is applied to the total sample, rather than to individual 
samples. This is the reason why the mean luminosities of 
the morphological samples are somewhat different from the 
observed values. We find that the M260/L ratios are inde- 
pendent of galaxy luminosity, morphological type and envi- 
ronment, consistent with the observational results. The pre- 
dicted amplitude of M260 /L is about two times the observed 
value, indicating that the model over-predicts the halo mass 
(or under-estimates the luminosity). 

As shown above, the predicted mass-to-light ratios in 
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Figure 10. The M260/L ratios for various samples taking into 
account the observational selections. The mass/light ratios are 
measured separately for spiral and elliptical galaxies, galaxies in 
low- and high-density regions, and all galaxies. The results of 
various simulation samples are connected with a solid line. For 
comparison, the results of M2002 are shown by symbols connected 
with a long dashed line. 



the simulations are about two times that of the observed 
value. Should we take this discrepancy seriously? In the GIF 
catalogs, model parameters are calibrated to reproduce the 
observed Tully-Fisher (TF) relation. This calibration de- 
pends on how the observed rotation velocities are related 
to the halo circular velocities and is hence uncertain. Kauff- 
mann et al (1999) have shown that their model underpredict 
the number density of L* galaxies by a factor of 2, and so a 
higher luminosity for a given halo mass may make the model 
prediction in better agreement with the observed luminosity 
function around L». However, assigning much higher lumi- 
nosities to dark haloes may lead to too high a luminosity 
density for the local universe (Kauffmann et al. 1999; see 
also Yang, Mo & van den Bosch 2002). 



4.2 Deriving halo masses from the GMCFs 

The relation between the aperture masses M260 and the 
mass of a galaxy's halo is complicated. We have already high- 
lighted one important problem- for "non-central" galaxies 
in groups, M260 will have almost no relation to the mass 
of the surrounding dark halo. Even if one were to impose 
an isolation criterion on the sample and select only central 
galaxies, the relation between M260 and M^aio is not trivial. 
For luminous central galaxies, M260 only measures the mass 
in the inner regions of their halos, while for faint central 
galaxies, M260 measures the mass beyond the virial radius 
of the halo . This is shown in Fig. [TIJ As one can see, M260 
is lower than the mean halo mass for bright central galax- 
ies and is higher than the mean halo mass for faint central 
galaxies. Over a large range of luminosity, the aperture mass 
M260 scales linearly with galaxy luminosity, while the halo 
mass scales as L . 

Fig. ^ shows the relation between Afhaio and M260 for 
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Figure 11. The aperture mass M26O (squares) and halo mass 
.Mhalo (circles) recovered from the GMCFs, together with the 
mean halo mass directly measured from the simulations (his- 
togram). Solid symbols are for CG samples and open symbols 
are for low-density samples taking into account of observational 
selection. The two straight lines show power laws M oc L 1,5 and 
M oc L, respectively. Note that over a large luminosity range, the 
aperture mass M26O in the simulation increases linearly with L, 
while the halo mass increases roughly as L 1,5 . 

central galaxies in the simulation. Once again we see that 
A/haio and M26Q are tightly correlated, but the relation is 
not linear. For haloes more massive than lO 13 ft" 1 M , M ha i 
is larger than M260, while for less massive haloes, Mhalo is 
smaller than h'heo- This relation can also be obtained using 
a simple analytical model. From numerical simulations, it is 
known that the concentration of a halo c is correlated with 
M. We use the c - M relation obtained by Bullock et al. 
(2002): 

/ M \ ~ ' 13 

c(M).9(-) , (6) 

where M„ is the characteristic non-linear mass scale at which 
the rms of the linear density field is equal to the critical 
overdensity for spherical collapse, a(M,) = S c ~ 1.69. The 
definition of the halo mass M is somewhat different from 
ours, but this relation changes only a little. With this model, 
we can predict Mhalo as a function of M260 ■ The solid curve 
in Fig. h2j shows this prediction. 

We now ask how one can derive halo masses from the 
full GMCF. One possibility is to use the results of the NFW 
fit to the GMCF to obtain an average 'halo mass' (denoted 
by Mhaio) for a set of lens galaxies that have no close neigh- 
bours of comparable brightness and are thus likely to be 
"central" objects. In Fig. [Tl] we plot the mass recovered 
this way as a function of galaxy luminosity for the central 
galaxies in our simulation (solid dots). A comparison with 
the mean halo mass measured directly from the simulation 
shows that the agreement is very good over a large range of 
halo masses. Thus, for central galaxies, this method works 
well in recovering the mean halo mass. The open circles in 
Fig. [n] show the recovered halo masses for the sample of 
galaxies in low-density environments defined in section 3.3. 



Figure 12. The MiaZo _ M260 relation. The solid line is the pre- 
diction of a simple model discussed in the text. The squares are 
inferred from the GMCF for CG samples. 

The results are very close to that for central galaxies, except 
at the very bright end where the halo mass is lower because 
bright galaxies in massive, rich clusters are missing. Thus, it 
is also possible to obtain the halo mass of galaxies by using 
the GMCFs of galaxies in low-density environments. In the 
next subsection we will apply this method to the data of 
M2002 to derive halo masses for the observed galaxies as a 
function of luminosity. 

The second method for deriving halo masses is to make 
use of the relation between Mhaio and M260 obtained from 
simulations (Fig. ^). This will be described in more detail 
in the next section. 

4.3 Halo masses of the observed galaxies 

We now apply our two methods of estimating halo masses 
to the actual observational data. Fig.[l^ shows the results 
of fitting a NFW profile to the GMCFs of M2002 (their 
Fig. 15). As in M2002, we have used only data points at 
projected radii smaller than 260/i _1 kpc. In order to obtain 
stable results, we have assumed that the concentration c 
depends on halo mass as in equation (^) . The open circles in 
Fig.jjr] show the halo mass obtained in this way as a function 
of galaxy luminosity. The errorbars on the derived masses 
assume that the observed errors on the GMCF at different 
radii are independent and are Gaussian distributed. Since 
the observed errors at different radii are actually correlated, 
the real errorbars on Mhaio should be smaller than those 
shown in the figure. 

The other way to derive Mhaio is to convert the aperture 
mass M260 into a halo mass using the relation shown in 
Fig.jl2|. These results are shown as solid squares in Fig.pr|, 
and the errorbars are obtained directly from the errorbars 
on the observed M^eo- The halo masses derived from these 
two methods are very similar. 

As one can see, the halo mass increases with luminosity 
roughly as a power law. A simple linear regression of all the 
data points gives 
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Figure 13. The NFW fits to the observed z' band GMCF data 
in Fig. 15 of M2002. Here we assume that the concentration c de- 
pends on the halo mass as Eq. (^) , and use T200 as the only free pa- 
rameter. The fits only extend over radii smaller than 260 /i _1 kpc, 
as in M2002. 
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Mean Luminosity (10 10 h~ a L o ) 



Figure 14. The derived halo masses M^ a i as a function of lu- 
minosity. The squares show M halo estimated using the relation 
shown in Figjl2[ The circles are the halo masses obtained from 
the fits in Fig. 16. Errorbars are estimated using 1-cr variances of 
the data points. To avoid confusion, the results of the NFW fit 
are shifted to the left. The solid line is a power law fit of the data 
points shown. 



2.o x io 12 /r 



*M f 



(7) 



where L, « 2.56 x 10 h~ Lq is the characteristic luminos- 
ity in the z'-band (Blanton et al. 2001). This scaling is sim- 
ilar to that obtained from the GIF simulation (see Fig. |ll|) . 
In terms of mass-to-light ratio, we can write 



M h 



75h 



(8) 



This mass-to-light ratio is about two times smaller than the 
simulation results. The origin of this discrepancy has already 
been discussed in Subsection 4.1. 

The power index 1.5 we obtain for the M^io-L relation 
is similar to that derived by Guzik & Seljak (2002) using 
an analytic model of the galaxy distribution in dark haloes. 
However, our value of the mass-to-light ratio for L» galaxies 
is about two times larger than theirs. One possible reason 
for this discrepancy is that Guzik and Seljak assume that 
galaxies in groups and clusters have the same halo masses 
as isolated galaxies of the same luminosity. If galaxy halos 
in groups are tidally truncated, their assumption will cause 
the halo masses of isolated galaxies to be underestimated. 
Current high-resolution N-body simulations (e.g. Klypin et 
al. 1999; Moore et al. 1999; Springel et al. 2001) show that 
subhalos of dark matter in galaxy systems have radii typi- 
cally much smaller than 200 kpc, and so their contribution 
to the lensing signals on the scales we are concerned here 
should not be important. 

Since the ratio M260/L does not depend on morpho- 
logical type, we may apply our method to spirals to check 
whether their halo masses are consistent with those de- 
rived from the Tully-Fisher relation (e.g. Hudson et al. 1998; 
Smith et al. 2001; Wilson et al. 2001; Seljak 2002). 

In the /-band, the observation of Giovanelli et al. (1997) 
gives 



Li 



1.0 x 10 1 



/ VpbB V 

U50 kms- 1 / 



h L 



ei ■ 



(9) 



where V a h B is an observational measure of the rotation ve- 
locity, usually taken as the maximum rotation velocity. The 
typical z' — I colour of a spiral galaxy is similar to that of 
the Sun, so its luminosity in solar units will have the same 
value in different bands. Thus, for an L» galaxy with a'-band 
luminosity 2.56 x lO lo /i -2 L0, the Tully-Fisher relation im- 
plies Vobs ~ 205kms _1 . If we include an internal extinction 
of about 0.3 m (see Verheijen 2001) this will increase V h s to 
about 225 kms" 1 . The halo mass for an L, galaxies obtained 
from the Tully-Fisher relation is thus 



M h 



4.8 x 10 12 h 



.1 rVo_ 
V V^bs 



M . 



(10) 



where V c is the circular velocity of the halo. We have as- 
sumed Slo = 0.3, and defined the halo radius such that the 
mean density within it is 200 times the mean density of the 
Universe. Comparing this with the lensing result given in 
equation (f?|), we see we require V hs/V c ~ 1.34. This ratio 
is consistent with the predictions of disk formation mod- 
els (e.g. Mo, Mao & White 1998; Mo & Mao 2000). The 
Vobs/Vc ratio obtained here is lower than that obtained by 
Seljak (2002) based on the mass determination published in 
Guzik & Seljak (2002). This is because our estimated halo 
mass is higher than theirs. 



5 DISCUSSION AND SUMMARY 

In this paper we have used galaxy catalogues constructed 
from a high-resolution numerical simulation to investigate 
the properties of the projected galaxy-mass correlation func- 
tion (GMCF) measured from weak galaxy-galaxy lensing 
observations. A comparison between the predicted GMCFs 
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with those obtained recently by McKay et al. (2002) from 
the Sloan Digital Sky Survey (SDSS) data shows that the 
observed dependence of the GMCF on galaxy luminosity, 
morphological type and environment is well reproduced in 
the model. In the simulations, the ratio between the aper- 
ture mass within a radius of 260fr~ 1 kpc and lens galaxy 
luminosity is independent of galaxy luminosity, type and 
environment. The value of M260/L is approximately equal 
to the ratio between the halo mass (defined to be the mass 
within the halo virial radius) and luminosity for L* galax- 
ies. M260 overestimates the halo mass for lens galaxies with 
L«X t and underestimates the halo mass for 

Using the simulations, we show that halo masses can be 
recovered by fitting the full GMCF with realistic halo den- 
sity profiles. If we apply this method to the observational 
data, we find that the halo mass scales with galaxy lumi- 
nosity as Mhaio oc L 1 ' 5 . The mean halo mass we derive for 
an L* spiral galaxy is consistent with that inferred from the 
observed Tully-Fisher relation. 

We find that the behavior of the GMCF can be signifi- 
cantly affected by satellite galaxies in high-density environ- 
ments. The observed GMCF can only be used to probe the 
properties of individual galaxy haloes if galaxy samples are 
limited to central galaxies. We have demonstrated that an 
effective way of doing this is to select bright and relatively 
isolated galaxies. In the future when lensing signal can also 
be detected for galaxies fainter than L, (many of which are 
expected to be satellite galaxies), detailed modelling of the 
galaxy distribution in dark halos will be required to compare 
models and observations. 

Our results suggest that the GMCF obtained from 
galaxy-galaxy lensing can provide important constraints on 
galaxy formation models. With the completion of the SDSS, 
such constraints will become an integral part of our under- 
standing of galaxy formation. 
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APPENDIX A: THE PROJECTION OF NFW 
PROFILE 

The properties of the projected NFW profile can be found 
in, e.g., Bartelmann (1996). The surface mass density can be 
expressed as a function of a dimensionless radius x = R/r s : 



E(aO = A f(x) 
where 

f(x<l) = 

f{x=\) = 

f(x > 1) 
and 



(Al) 



In 



(A2) 



x - 1 



1 - 



atanVs 



Vx 2 ~ 1 



A = 2 r s 8 p c . 

The mean surface mass density within radius x is: 

E(< x) = A g{x) 

where 



(A3) 
(A4) 



2 I, I In 

9(x<l) = -5 In - + 



g(x = l) = 2 + 21n : 



(A5) 



, , . 2 / , x atanv^ 2 — 1 

9{X>1) = ?h + ^=r 

The mass density contrast is simply 
AE+(x) = A[g(x) - f{x)]. 



(A6) 



APPENDIX B: CONTRIBUTION BY 
SATELLITE GALAXIES 

If the lens galaxy is not located at the center of its halo, the 
average surface mass density and the GMCF will depend on 
an extra parameter R c , the distance between the galaxy and 
the center of the halo. The average surface mass density in 
this case is 



(Bl) 



E ' (r) = ^/ E{r ' )de 

where S(r') is the NFW surface mass density and r' = 
VR'i + r 2 + 2R c rcos6. 

In Fig. Bl we show the surface mass density (upper 
panel) and the density contrast (lower panel) as a function 
of radius R for various choices of R c . The parameters are set 
to be r2oo = 1.2/i _1 Mpc, c = 6.0, and R c changes from 
to 0.4 h~ Mpc. The figure shows that the surface mass den- 
sity has a peak near R c , while the density contrast reaches a 
minimum near R c and then increases rapidly with increas- 
ing radius. Clearly, the behavior of the GMCF for satellite 
galaxy is very different from that for central galaxy (corre- 
sponding to R c — 0). The complex behavior of the GMCF 
of satellite galaxies complicates the interpretation of the ob- 
served GMCF, and so it is important to separate central 
and satellite galaxies in the galaxy-galaxy lensing study. As 
we have demonstrated above, selections of central galaxies 
can be made quite effectively by using bright and relatively 
isolated galaxies. 
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